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Abstract - When ignorance due to the lack of knowledge, 
modeled as epistemic uncertainty using Dempster-Shafer 
structures on closed intervals, is present in the model pa- 
rameters, a new uncertainty propagation method is neces- 
sary to propagate both aleatory and epistemic uncertainty. 
The new framework proposed here, combines both epistemic 
and aleatory uncertainty into a second-order uncertainty 
representation which is propagated through a dynamic sys- 
tem driven by white noise. First, a finite parametrization 
is chosen to model the aleatory uncertainty by choosing 
a representative approximation to the probability density 
function conditioned on epistemic variables. The epistemic 
uncertainty is then propagated through the moment evolu- 
tion equations of the conditional probability density func- 
tion. This way we are able to model the ignorance when the 
knowledge about the system is incomplete. The output of the 
system is a Dempster-Shafer structure on sets of cumulative 
distributions which can be combined using different rules of 
combination and eventually transformed into a singleton cu- 
mulative distribution function using Smets' pignistic trans- 
formation when decision making is needed. 

Keywords: Uncertainty Propagation, Epistemic Uncer- 
tainty, Aleatory Uncertainty, Dempster-Shafer, Moment 
Evolution, Pignistic Transformation, Ignorance 

1 Introduction 

The paper presents a novel framework to characterize the 
response of stochastic dynamic systems with model param- 
eters governed by epistemic uncertainty. We distinguish be- 
tween aleatory uncertainty which arises due to the stochas- 
tic behavior of the system and epistemic uncertainty which 
is used to model the ignorance in the model parameters. In 
this last case we do not know if the parameter in discussion 
is random or not, and if it is random what is its underlying 
probability distribution 1161 . 

Ferson and Ginzburg ||5| argues that when dealing 
with propagating uncertainty through mathematical mod- 
els, different calculation methods are required to propagate 
aleatory and epistemic uncertainty. However, the segrega- 
tion of both types of uncertainty during propagation, as well 



as the accurate propagation of either one of them are two 
challenges in creating a method to propagate both epistemic 
and aleatory uncertainty. 

Finding only the response of a stochastic dynamic sys- 
tem with no uncertainty in the parameters, is a very ac- 
tive area of research. The problem in this case is to solve 
for the evolution of the probability density function (pdf), 
pit, x), corresponding to the state x of the dynamic system. 
The evolution of the pdf is governed by the Fokker-Planck- 
Kolmogorov (FPK) equation with analytical solution known 
only for stationary pdfs of a limited class of dynamic sys- 
tems Q [TSl . Different approximate methods exists in the 
literature to represent the pdf with a finite number of pa- 
rameters such as finite-difference techniques, Monte Carlo 
methods, Gaussian Mixtures models 1201 and Gaussian clo- 
sure methods. 

These methods provide approximations to the pdf when 
the prior pdf is precisely known, the process noise is per- 
fectly characterized and the parameters in the model are ex- 
actly known. This is also known as the Bayesian dogma of 
precision ll22l . however, in practice, these precise values are 
difficult to obtain due to the amount of information avail- 
able, incomplete knowledge of the system or systematic un- 
derestimation of uncertainty which arises in the elicitation 
process. This becomes of great importance when doing haz- 
ard risk assessment and decision making. 

The framework proposed here, propagates both epistemic 
and aleatory uncertainty through a dynamic system driven 
by white noise using a two-level hierarchical model. On 
the first level we model the aleatory uncertainty by choos- 
ing a representative approximation for the pdf with a finite 
number of parameters. In this paper we model the pdf con- 
ditioned on epistemic variables as a Gaussian distribution. 

On the second level is the epistemic uncertainty that we 
have in the moments of the stochastic solution, in this case 
the first two moments. The epistemic uncertainty is modeled 
in this work using finite Dempster-Shafer structures where 
the focal elements are closed intervals. The closed intervals 
are propagated through the moment evolution equations us- 
ing a recent proposed method based on Polynomial Chaos 



- Bernstein Form 1211 . Thus the framework provides the 
means to propagate both epistemic and aleatory uncertainty 
separately Q as well as combine them when a decision has 
to be made. 

After the propagation we obtain finite Dempster-Shafer 
structures, only this time the focal elements are sets of dis- 
tributions represented by probability boxes [61 . The struc- 
tures obtained characterize the uncertainty in the quantity of 
interest, namely the cumulative distribution function (cdf) 
of the response of the system, and they be used with classic 
decision theory by using the pignistic transformation [191 
whenever a decision is required. 

In Section|2]the problem to be solved is defined, followed 
by a background of Demspter-Shafer structures on closed 
intervals and decision making under ignorance in Section 
[3] The approximate uncertainty propagation framework is 
developed in Section |4] with resolutions for the two posed 
problems. In Section |5] a proof of concept example is pre- 
sented and the conclusions and future work are discussed in 
Section |6] 

2 Problem Statement 

Consider the following first order model with stochas- 
tic forcing and uncertain initial condition modeled using a 
Gaussian pdf: 



i(t) 



f{x,oi,t) 



A/'(xo; ^0,0-0 ) 



(1) 
(2) 



The Gaussian white noise process r(i) has the autocorre- 
lation function £'[r(t)r(T)] = q^S{t — t), and a is a vector 
of 71 parameters. 

The variables of the model can be segregated in two 
sets: the set of aleatory variables given by a{t) = 
{x{t),r{t)} and the set of independent epistemic vari- 
ables given by e(i) = [es,eD{t)], where e£)(i) = 
{TOi(i), 7712 (i) . . . TOfc(t)} is a time-variant vector of epis- 
temic moments used to characterize the probability distribu- 
tion of the response, p{t, x\eD), and e^ = [a, q^] is a time- 
invariant vector of epistemic variables. Given the initial con- 
dition, the first two moments are given by ?Tii(0) = po and 

?Tl2(0) =(J^+ pI- 

If a part of the model parameters, a, are characterized by 
aleatory uncertainty then they can be treated as additional 
state variables and they can be augmented to the aleatory 
vector a(t). The evidence about the model parameters in 
e(t) is modeled here using Dempster-Shafer structures on 
closed intervals and the uncertainty of the aleatory variables 
in a(i) is quantified by probability distributions. 

When the model parameters are known, the response of 
the system is characterized only by aleatory uncertainty 
which is represented here by the cumulative distribution 
function. 



where eg = e(0) = [es,e£)(0)] and p(i,.T|eo) is obtained 
in general by solving the FPK equation. 

Given the uncertain model parameters, we are interested 
in solving the following three problems: 

1. The main objective is to find the induced Dempster- 
Shafer structure about the response of the system, 
■m[P{X < Xf\eo)]. Hence, we are looking to deter- 
mine the focal elements of the response as well as their 
corresponding probability masses. 

2. Construct a cumulative density function using the pig- 
nistic transformation in order to use the expected utihty 
theory for decision making. 

The last problem represents a secondary objective and is 
included here in order to present the practicality of the ap- 
proach. 

3 Background 

A DS-structure on closed intervals is a collection of 
interval-valued focal elements and their associated basic 
probability assignments 161 . We write the uncertainty of 
the epistemic variable x is represented by the following DS- 
structure as: 

X^ { {[Xi,Xi],pi) , {[X2,X2],P2) , ■•■ {[x„,Xn],Pn) } (4) 

Given two nondecreasing functions F and F_, where 
F, F : R -> [0, 1] and F{x) < F{x) for all x G R, we 
can represent the imprecision in the cumulative distribution 
function, F{x) = P(X < x), by the probability box (p- 
box) [F,F] as follows: F{x) < F{x) < F{x) Jgl. 

A Dempster-Shafer structure on closed intervals can in- 
duce a unique p-box, while the inverse is not uniquely deter- 
mined. Many Dempster-Shafer structures exist for the same 
p-box. Given the following body of evidence, Eq.©, the 
cumulative belief function (CBF) and the cumulative plau- 
sibility function (CPF) are defined by: 



CBF{x) = F{x) = Y. 



Pi 



X.i<X 



CPF{x) = F{x) = ^ P, 



(5) 
(6) 



Thus the cumulative distribution function is bounded as 
follows: 



CBF{x) < P{X <x)< CPF{x) 



(7) 



In order to compute expectations using the interval-based 
belief functions, one needs to build a probability density 
function given the belief structure. Based on Smets' pig- 
nistic transformation one can define the pignistic probability 
density function pBet {x) WT\ as a finite mixture of continu- 
ous uniform distributions: 



P{X<Xf\e^) 



^f 



p{t, x\eo)dx 



(3) 



PBet 



(-)=E: 



Pi 



-I{x,[Xi,Xi]) 



(8) 



where the indicator function I{x, [xj, Xi]) is given by: 



I{x,U,^x,]) = 



1, for X e [x^,Xi] 



0, 



otherwise 



Using the above defined pignistic probabihty density 
function pBet{x) one can compute the expectations needed 
in the decision making process. Example the expected value 
of X given the pignistic probability density function si given 
by: 

/ + 00 1 '"■ 

XpBet{x)dx = - ^Pife +Xi) (9) 

The pignistic transformation constructs a singleton pdf 
and makes the expected utility theory applicable, however 
it ignores the ignorance ||9l- One possible way to incorpo- 
rate the ignorance into the decision process is to construct 
a scalar measures which quantifies the total amount of ig- 
norance and take a decision only if the level of ignorance is 
relatively low. Define the normalized integral of the degree 
of ignorance: 

1 r 

NIDI = / [CPF{x) - CBF{x)]dx (10) 

where x = iiiini<i<„a;j and x = maxi<i<„a;,;. The 
NIDI is a scalar measure with values ranging between 0.0 
and 1.0 and it summarizes the confidence in the pignistic 
pdf. A value equal with 1.0 denotes that we are dealing with 
interval uncertainty where we only know the bounds of the 
variable, and a value closed to 0.0 means that we know the 
pdf of X precisely. Thus, one can be more confident in ap- 
plying the expected utility theory if the amount of ignorance 
if low (eq. less than 0.1) and can try to defer the decision 
making and gather more evidence if the level of ignorance 
is not tolerable. 

4 Proposed Approach 
4.1 Time evolution of moments 

In order to keep both types of uncertainty segregated dur- 
ing propagation, one can use the moment evolution equa- 
tions associated with Eq.©, which is a set of ordinary dif- 
ferential equations (ODE) governed only by epistemic un- 
certainty. Thus, after the applying the Ito's lemma to Eq.© 
ifTTl . the time evolution of the moments conditioned on the 
epistemic variables is given by 

'^'^^^^^^'"'"^ = -n^.{x)f{x,c.,t)\e^] + iE[5V.4^)|eo] (11 

where ip{x) = x^ . 

The nonlinearity of the model is assumed to be polyno- 
mial, however the general case can be handled using numer- 
ical integration schemes such as Gaussian quadrature. Thus, 
when the nonhnearity in the model is polynomial, 



the final form for the evolution of the moments conditioned 
on the epistemic variables, is given by 

" 1 

mfc|eo = -fc^aimi+fc_i|eo + -k{k ~ l)g^mfe_2|eo . (13) 

1=1 

Observe, that in Eq.(fT3l) the lower order moments depend 
on the higher order moments, thus we have an infinite hier- 
archy of moment equations. To truncate the infinite chain, 
one can use moment closure schemes which assume a spe- 
cific class of probability distributions to express higher order 
moments in terms of lower order ones. 

In this paper we are using the Gaussian closure method, 
hence we are interested in propagating only the first two mo- 
ments, rrii and m2, under the assumption that p(<, a;|e£)) ~ 
N{x] /z, (T^|e£)). To express the moments of order fc > 2 in 
terms of the lower order moments one can use the following 
relation [T]: 

mfe = E ((X - ^if) - Y,{-lf~' r\ m,/-' (14) 

where mo = 1, ji = mi, a^ ~ m-i — mf, and the central 
moments of the Gaussian distribution are given by 



E((X-/.f) 



/c! 7 • 

2'=/2(fc/2)! ' ^^ even 
0, k is odd 



After using the Gaussian closure, the time-variant vec- 
tor of epistemic variables is given only by the first two mo- 
ments, eu{t) ~ [mi(i),TO2(i)]^, and the system of ODEs 
( fT3] l is transformed into the following system: 

eo = g(eD,es) (15) 

eD(0)~{([eoi(0),ei,i(0)],pDi), 

([£02(0), 602(0)], PC2), ... } 

es ~ { ([esi,esi],psi) , ([§52, es2],PS2) , . . . } 

Under the Gaussian closure assumption, p{t,x\eQ) = 
p{t, x\e£)) where e^i is obtain given Gq and the Eq. ( fTSl l. 

We are interested in finding the induced Dempster-Shafer 
structure in e£)(t) at time t > 0. Given the knowledge 
about elements of Gq as independent DS structures on closed 
intervals, we need to find the belief structure associated 
with e£)(t). The DS-structure describing the uncertainty for 
eoit) can be obtained using Yager's convolution rule for 
DS-structures under the assumption of independence 1241 : 



m[eDit)]{C,j) 



E 

eD=g(ei3,es) 

ei,(0)eAi 

esSB, 



m[ei5(0)](A,)m[e5](Bj) 



PS] 



f{x,a,t) = y^ ajx"- 



(12) 



(16) 
where Cij = [eo^j{t),eDijit)],Ai = [e^.^{0),emiO)] and 
Bj = [eg , esj]. Here the only unknown is the collection 
of focal elements Cij ofeoit) . 

Thus the problem of finding the mapping of a body of ev- 
idence on closed intervals is reduced to interval propagation 
ifTSl . This problem can be solved using the advanced tech- 
niques developed in the interval analysis field ITOl . How- 
ever, due to the dependence problem the obtained bounds 



are conservative which is detrimental to the belief structure, 
since the evidence is assigned automatically to other ele- 
ments which are not in the body of evidence. This problem 
becomes more acute when the uncertainty has to be propa- 
gated over a period of time. 

The dependence problem in interval arithmetics can be 
avoided by using Taylor models and the remainder differen- 
tial algebra for bounding the range of the response of ordi- 
nary differential equations under both initial value and para- 
metric interval uncertainty ifTSl l2l [T3JI . In Ref.H are pre- 
sented several examples in propagating uncertainties repre- 
sented by probability boxes using Taylor model methods and 
interval arithmetics. Also an excellent review on interval 
methods for initial value problems is presented in Ref . |[T4| . 

A method based on polynomial chaos and the Bernstein 
form is presented in Ref. [21] to propagate DS-structures on 
closed intervals through nonlinear functions. Given a func- 
tion of random variables with compact support probability 
distributions, the intuition is to quantify the uncertainty in 
the response using polynomial chaos expansion and discard 
all the information provided about the randomness of the 
output and extract only the bounds of its compact support. 
To solve for the bounding range of polynomials, we have 
proposed to transform the polynomial chaos expansion into 
the Bernstein form, and use the range enclosure property of 
Bernstein polynomials to find the minimum and maximum 
value of the response lO- 

The PCE is mathematically attractive due to the func- 
tional representations of the stochastic variables. It separates 
the deterministic part in the polynomial coefficients and the 
stochastic part in the orthogonal polynomial basis. This be- 
comes particularly useful in characterizing the uncertainty 
of the response of a dynamical system represented by ordi- 
nary differential equations with uncertain parameters such 
as in Eq.(fT5ll. The result is a set of deterministic differen- 
tial equations which can be solved numerically to obtain the 
evolution of the polynomial coefficients. 

4.2 Interval Uncertainty propagation 
through ODEs using PCE and the 
Bernstein Form 

Given the initial value problem in Eq.(fT6]l with interval- 
valued initial condition and parameters. 



Afc 



gkieD,es) 
eD(0)G[A,A] 
esG [B,B] 



(17) 



our goal is to approximate the range of the state variables e^ 
at time t, where e^ is the fc-th component of eD{t). Here, 
intervals of general type have been used in order to reserve 
the indexes for later use. 

The problem can be transformed into finding the stochas- 
tic response under the assumption that both the initial con- 
dition and the model parameters are uniformly distributed. 

Thus we define e%{0) ~ U{A\A') and eJg - U{B;\W), 



where e^ (0) and e^. is the i-th component of en (0) and the 
j-th component of e^ respectively. 

We expand both the uncertain arguments and the response 
of the system using the finite dimensional Wiener-Askey 
polynomial chaos ||23ll : 

p-i 

4 = E4^j(^') ^h^^'^ 6 ^^(-1,1) (18) 

p-i 



^D = T.^D^M^) Where p^t±P}l 



i=0 



r\p\ 



(19) 



Here r is the number of uncertain input variables and is 
equal with the sum of the size of the vector 6^3(0) and e^, 
and p is the order of the polynomial expansion. The basis 
function ipj is the j-th degree Legendre polynomial and ■0,; 
is a multidimensional Legendre polynomial and the poly- 
nomial coefficients are initialize such that they match their 
initial uniform assumption. 

We are interested in finding the polynomial coefficients 
e|jj of the response after t sec. Substituting Eqs.(fr8]l-(fT9]l in 
Eq.dTTJi and using the Galerkin projection and the orthogo- 
nality property of the polynomials one obtains a system of 
P deterministic differential equations which can be solved 
numerically to obtain the PC expansion coefficients of the 
fc-th moment. 



-Di 



< gk,ipk > 



(20) 



where < . , . > represents the inner product operator and 
can be evaluated in general using sampling or quadrature 
techniques. In this particular case since the nonlinearity is of 
polynomial type, the inner product between different Legen- 
dre basis functions can be computed a priori, speeding this 
way the numerical integration of the ODE. 

After t sec, by integrating Eq.(|20li. on obtains the poly- 
nomial coefficients which define the stochastic response of 
the system. However, we are only interested in finding the 
bounds that enclose the response of the system. In Ref. 1211 
it is shown that by bringing the polynomial chaos expansion, 
Eq.(fT9ll, to a Bernstein form using the Garloff 's method [S], 
one can efficiently find the range of the compact support 
thanks to the enclosing property of Bernstein polynomials. 

4.3 DS structures on Probability Sets 

By propagating the DS structures through the moment 
evolution equations we obtain, at time t, an induced DS 
structure for the moments that characterize the probability 
distribution of the response. Since in this paper we have cho- 
sen to approximate the conditional pdf of the response using 
Gaussian density functions, we are only interested in the DS 
structures for the first two central moments: the mean, /^, 
and the variance, a^. Thus, using Eq.dTSl) the following DS 
structure is obtained for the first two moments: 

^D{t) - { ([e£,i(i),eDi(t)],p£,i) , 

{[BD^{t),eD2{t)],PD2) , ... } (21) 



For each focal element we obtain a pair of two intervals 
that bound the range of the mean and the variance of a Gaus- 
sian density function. 

Consider now the following Gaussian density func- 
tion, J\f{x; fi,a'^), with uncertain parameters, fi and a^, 
given by two intervals: [fJ.,]l] and [af^, cr^]. Let us de- 
note the cumulative distribution function N {x f, jj,, a'^) = 
jy^Af{x;ii,a^)dx. In Ref.|125j| it is shown that all the 



normal cdfs are bounded by two functions, N.{xf) < 
N{xf; ^, CT^) < N{xf), that can be computed analytically. 



N{xf) = 



Njxf, u,a ) Xf > jj. 

N{xf;fj.,a'^) Xf<^j. 

N{xf;Jl,a^) Xf >J1 

N{xf\JL,g^ Xf <Jl 



(22) 



Using the above envelope property and the DS structure 
in Eq. (l2n i. we obtain the following induced DS structure for 
the response of the system. 



P{X < Xf\eo) 



\N^{xf),N,{xf)],pni), 

{[N^{^f),N2{xf)],PD2) , 



(23) 



Thus, we are modeling a system with second order uncer- 
tainty and our credal set is defined as a DS structure over 
p-boxes. Having a structure like this makes the problem of 
decision making difficult, since we are dealing with many 
envelopes of cdfs instead of a singleton cdf or just an im- 
precise probability represented using a p-box. However, the 



following subsection presents how the pignistic transforma- 
tion can be used whenever decision making is needed. 

4.4 Constructing a singleton CDF and Deci- 
sion Making using DS structures on Prob- 
ability Sets 

In order to compute expectations using the interval-based 
belief functions, one needs to build a probability density 
function given the belief structure. Based on Smets' pig- 
nistic transformation one can define the pignistic probability 
density fVinc\ion pBet{x) ifTTl as a finite mixture of contin- 
uous uniform distributions. Thus, given the DS structure in 
Eq.(l23]l. we can construct a singleton cdf using the pignistic 
transformation as shown in Eql24l 



PBet{X<Xf) ^ EBet\P{X<Xf\eo)\ (24) 

1 "" / — \ 



4 = 1 ^ ^ 

Here, n^) is the total number of focal elements in the body 
of evidence Eq.(l23]). Notice that this is different from the 
law of total probability where a probability density function 
is constructed for the epistemic variables either eo or bd- 
In both cases the pignistic transformation is applied much 
earlier in computing the quantity of interest. 

Furthermore one can construct an ignorance function by 
applying Eq.(fTOb for each xj. This associates a point-wise 
measure of confidence in constructing the pignistic cdf. 



IgF{xf) 
where 



N(xf) 



N{xf) 



CPFp(x<xf\eo){z) — CBFp(x<x;\eQ){z) 



Az 



(25) 



iV(2^-/) = , min Kiixf) and iV(a;/) = max N,{xf) 

CPFp(X<X;\eo){z) = 5Z PO* 

iV;{a:j:)<z 

CBFp^x<xi\Bo){z) = ^ Voi 

7fi{xf)<z 



Similar to Eq.dTOli. we can construct a scalar measure to 
summarize the total amount of ignorance by integrating the 
above ignorance function. 



NligF -■ 



1 



IgF{xf)dxf 



(26) 



Given the pignistic cdf, a;,ni„ and aJmax are the pth and 
(100 — p)th percentiles respectively (eq. p ^ 0.05). The 
NligF is a number between 0.0 and 1.0 and can be used 
similarly as NIDI, in Section [3] to make decisions using 
the pignistic cdf in Eq. (|24] |. if it's value is small (eq. less 
than 0.1) otherwise defer the decision if it is permitted and 
gather more evidence. Again, a value of 1.0 denotes that 



we are deaUng with interval uncertainty and a value of 0.0 
means that we now the cdf precisely. 

5 Numerical Simulation 

Consider the following linear dynamic system driven by 
Gaussian white noise: 

x + aix = ri(t) (27) 

p{to,xo) = A/'(a;o; 1.1, 2.42) (28) 

where the autocorrelation function of the noise is 

£'[ri(i)ri(T)] = qfS{t — r), and both oi, qi are epistemic 



variable, described by the following body of evidence: 

ai ~ I ([0.86, 0.9], 0.2) , ([0.89, 0.96], 0.8) I (29) 
qi - I ([0.2, 0.3], 0.6) , ([0.3, 0.4], 0.4) I (30) 

We are interested in finding the quantity of interest P{t = 
2,Y < —0.5) as well as construct the entire cdf, P{t = 
2, r < y), after t = 2 sec. 

Here, the exact conditional probability density function is 
Gaussian, 



p{t, x\ai,qi) = 7V(x; fi{t), a'^{t)\ai,qi) 



(31) 



because the model is linear and the initial condition and the 
process noise are normally distributed. 

Thus, we are interested in finding the induced DS struc- 
tures for the moments of the normal distribution in Eq.(|3T|i. 
The time evolution of moments through which the epistemic 
uncertainty is propagated is given by: 



nil 

1712 



— 2aim2 

mi(0) = 1.1 and m2(0) = 2.42 



(32) 



^Z? 



Using the PCE - Bernstein form described in Section 
14.21 to propagate interval uncertainties through the ODE in 
Eq.(l32b. and solving for the central moments we find the DS 
structure for /i and a^ after i = 2 sec, shown here in Table 
[Hand graphically represented in Fig[l^. 



The DS structure over p-boxes is obtained by using the 
envelope property in Eq.(l22]i. The focal elements and their 
bpa for P{t = 2,Y < y\ai,qi) are presented in Fig[T}5. 
From this structure it is easy to obtain the DS structure for 
P{t = 2,Y < — O.Sjfli, qi) as it is shown in FiglT];. Now, 
using the pignistic pdf for P{t ^ 2,Y < — 0.5|ai, qi) one 
can construct an estimate for the quantity of interest using 
Eq.(l24b as well as provide a measure of confidence for this 
estimate: 

Psetit = 2, y < -0.5) = 1.07% with NIDI = 1.56% (33) 

Eq.(l24ll can be used to construct a singleton cdf for y to 
take an action using the expected utility theory. Also using 
Eq.jZSl) we can obtain an ignorance function to indicate the 
point-wise confidence in the constructed cdf as well as pro- 
vide a scalar measure of trust, Eq.(l26b, in using this cdf in 
decision making. Both the constructed cdf and the ignorance 
function are presented in Fig[T]F. 

Just for comparison purposes the probability theory has 
be used to obtain the total probability of y under the as- 
sumption that the epistemic variables vary randomly. Us- 
ing Laplace's Principle of Insufficient Reasoning and we can 
transform the epistemic variables ai and qi into the aleatory 
variables al and ql . For this, we use the pignistic transfor- 
mation to find the following pignistic pdfs: 

PBetK) = 0.2^^(0.86, 0.9) + 0.8Z^(0.890.96) (34) 

PBetiql) = 0.6^^(0.2, 0.3) + 0.4^(0.3, 0.4) 
With this transformation the uncertainty in the model in 
Eq.dZTJ) is represented solely by probabilities. Thus an esti- 
mate for the quantity of interest is calculated as in Eg [35] 



PMcit = 2,Y <~0.5) 



P(t = 2, y < -0.5\al , qi)pBet{al)pBet{qi)daldqi 



(35) 



/-0.5 
Af{x;fj,,a \al,ql)dx 
-oo 



= 0.85% 



PBet{al)pBet{qi)daldql 



In Fig[T^ it is plotted the histogram for P{Y < 
— 0.5[a*, ql) given 1 million samples drawn from the pdfs 
inEq.®. Also the empirical cdf for P(y < -0.5\al,ql) 
is shown in Fig[TJl and the total cumulative probability for y 
is presented in Fig[TJ. 

The probabilistic estimate for the quantity of interest dif- 
fers from the one given by the present approach due to the 
assumption that the model parameters vary randomly. Since 
we have used the pignistic transformation to obtain a pdf for 
the model parameters and then propagate this pdf through 
the dynamic system, the concentration in the probability 
mass seen in Fig[T^ cannot be explained given the evidence 
we have started with. This argument is nicely explained for a 
simple example by Ferson in Ref.|5j. More over the NIDI 
measure indicates how much can we trust our estimate given 



our lack of knowledge about the model parameters. The 
threshold that indicates if the estimate can be trusted can 
only depend on the magnitude of the consequences is the 
state of the system drops below —0.5. 

The two cdf for y in Fig[TJ are slightly different, however 
for the cdf constructed using the present method we also in- 
clude an ignorance function as a point- wise measure of con- 
fidence and a number, NligF, which similar with NIDI. 
Here, while the probabilistic approach offers no other alter- 
native just to use the expected utility theory to take an action, 
our approach through the use of the ignorance function and 
the NligF measure can postpone the decision making if the 
amount of ignorance is high. The use of ignorance function 
and NligF remains to be studied in the decision making 
problem, and is not the purpose of the current paper. 



Table 1: The induced DS structure for {/i, cr^} 



m„i ([0.86, 0.9]) =0.2 



mai ([0.89, 0.96]) =0.8 



1,1 ([0.2, 0.3]) =0.6 m{^,^2}n0.182, 0.197], [0.055,0.090]) =0.12 m{^,„2} ([0.161, 0.186], [0.047, 0.084] j =0.48 
1,1 ([0.3, 0.4]) =0.4 m{^,^2}M0.182, 0.197], [0.082,0.129]) =0.08 m{^,^2} ([0.161, 0.186], [0.076, 0.122] ) =0.32 
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Figure 1 : Numerical simulation: DS structures, pignistic cdfs and ignorance function after i = 2 sec 



6 Conclusions 

A new framework has been presented to propagate both 
aleatory and epistemic uncertainty through dynamic systems 
with stochastic forcing. This has been achieved by using a 
second-order uncertainty model to propagate both types of 
uncertainty. 

First an approximation to the probability density function 
of the response is assumed and it is used in the moment 
closure scheme to find the corresponding time evolution of 
the moments. The epistemic uncertainty is mapped through 
these equations and the final response of the system is mod- 
eled as a Dempster Shafer structure on probability boxes. 

The paper incorporates a previous work of the author to 
propagate interval uncertainties through ordinary differen- 
tial equations using the polynomial chaos expansion and the 
Bernstein form. However, other methods such as Taylor 
models and interval arithmetic can be used to propagated 
the focal elements through the moment equations. 

It is shown that the Demster Shafer structures obtained 
can be used in constructing estimates of quantity of inter- 
ests, such as probabilities of failure, and also in constructing 
entire cumulative density functions which can be used in 
decision making. In addition, an ignorance function and 
a scalar measure can be computed in order to hint the 
confidence that one can have in the estimates given the 
initial body of evidence concerning the model parameters. 
For now one may choose not to pursue with the decision 
making if the scalar measure is over a certain threshold. 
The use of the ignorance function in conjunction with the 
expected utility theory is set as future work. 
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